Quasiparticle transport in the vortex state of d-wave superconductors. 
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, We calculate the magnetic field and temperature dependence of the Raman response, superfluid 

■ density and the NMR relaxation rate in the vortex state of a d-wave superconductor arising from the 

Doppler energy shift of extended quasiparticle states. Our results are valid both at low temperatures, 
where we observe scaling with variable Tif -1 / 2 and obtain explicit form of the scaling functions, 
and beyond this region. We derive a universal frequency dependent scaling relation for the Raman 
response, and discuss the breakdown of the single relaxation rate approach to NMR response. 
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In the last few years there has emerged a consensus regarding the d-wave symmetry of the order parameter in the 
hole doped high-T c cuprates. While in conventional, s-wave, superconductors a finite energy gap exists everywhere 
^ on the Fermi surface, the d-wave order parameter has lines of nodes, which leads to a gapless excitation spectrum 
along certain directions in momentum space. Consequently the low temperature behavior of the thermodynamic 
9^, and transport coefficients in the high-T c materials are qualitatively different from those of the s-wave compounds. 
^5 ' Properties of the vortex state in "unconventional" superconductors also differ significantly from those obtained within 
the framework of s-wave superconductivity, and a clear understanding of these properties is essential for interpretation 
of the experimental results and a better grasp of the novel physics associated with the superconducting state of the 
high-T c materials. 

Volovik pointed out that while in an s-wave superconductor the density of states is determined at low fields 
by the quasiparticles in the vortex core, both the density of states and the entropy in d-wave superconductors are 
Q ■ dominated by the extended quasiparticle states which exist even at zero temperature in the nodal directions of the 
O , order parameter. A remarkable consequence of this behavior is that the specific heat of such a superconductor varies 
as VH rather than linearly in the applied field. Other authors || used the Dirac form of the low-energy excitation 
t-H , spectrum of nodal quasiparticles to demonstrate that thermal and transport coefficients exhibit scaling with TH^ 1 / 2 . 
However, strictly speaking, the analysis applies only to clean superconductors, and the energy spectrum is only Dirac- 
like at energies small compared to the gap amplitude Ao; a crossover to a different scaling regime followed by the 
breakdown of scaling have been predicted at (T /T c )(H c2 / H) 1 / 2 ~ 1 §]|. 

Very recently Kubert and Hirschfeld (JJ^] placed these scaling arguments in the framework of the Green's function 
formalism capable of treating both the energies of order of the gap and the effects of disorder. These authors argued 
that for a short coherence length superconductor the typical spacing of the energy levels in the core, A^/Ef, where Et 
is the Fermi energy, is large so that only one or a few states (if any) exist there, and suggested that the contribution 
of the vortex cores to the transport coefficients is negligible over a wide region of H and T. They proposed to 
account for the effect of the magnetic field on the extended states semiclassically by introducing a Doppler shift due 
to circulating supercurrents, which for H <C H c2 are approximated by the superfluid velocity field around a single 
i-^J ■ vortex v s = fiOjlmr, where r is the distance from the center of the vortex and 9 is the azimuthal angle in real 
space. The authors of Ref. |J investigated in detail the breakdown of scaling of both specific heat and the thermal 
conductivity with increased impurity scattering, and the results agreed remarkably well with recent experiments || . 

In this work we use the same approach to examine the effect of magnetic fields H c \ < H <C H c2 on the Raman 
response, superfluid density, Knight shift, and NMR relaxation rate in a d-wave superconductor. These are issues of 
considerable experimental interest: recently Blumberg et al. analyzed the changes induced by a magnetic field in 
the electronic part of the Raman response; NMR relaxation rates and the Knight shift are measured in fields of up 
to 10T j8j, in the vortex state, while muon spin rotation (/iSR) is used to determine the low temperature penetration 
depth in a magnetic field which is related to the superfluid density. In our analysis we obtain scaling similar 
to that suggested in Ref. H0| and give the explicit form of the scaling functions at low temperatures, however our 
results remain valid in a wider parameter range. We do not include impurity scattering in the present calculations, if 
included, it affects only the extreme low frequency (low temperature) part of the Raman (NMR) response, and the 
results presented here remain valid beyond this narrow region. 

We employ the single particle Green's function which is obtained by introducing the Doppler shift into the BCS 
function ju||| 
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where w n is the fermionic Matsubara frequency, (k is the energy of a quasiparticle with momentum k, measured with 
respect to the Fermi level, and Tj are Pauli matrices. The Green's function depends on the coordinate r in real space 
via the superfluid velocity v s . Now thermal and transport coefficients can be calculated using the standard approach 
jyj, however, they become local quantities which have to be averaged over a unit cell of the vortex lattice. Following 
Rcf. || we approximate this unit cell by a circle of radius i?, where 2R = ^ Q ^/2~na~ 1 {H C 2/ H) 1 / 2 is the intervortex 
spacing, and a is a geometric constant of order unity, so that the average of a quantity /(r) is given by 

Here we consider a model cylindrical Fermi surface and an experimental arrangement with the magnetic field parallel 
to the axis of the cylinder, Hj|c. 

Raman response. The Raman intensity is proportional to the imaginary part of the zero momentum density-density 
correlation function 

X (in n ) = - 7 |T ftfe)Tr [r 3 G(k, ^ m )r 3 G(k, iu m - iO n )] , (3) 

where 7_R/ S (k) is the Raman vertex. In Eq.(||) vertex corrections due to Coulomb screening have been ignored, they 
appear only in the fully symmetric channel (A\ g ), while we are interested primarily in the B\ g and Big scattering 
geometries. After doing the Matsubara sum and analytically continuing the response function to real frequencies, we 
obtain the local Raman response 

X "(n; r) = y R N(0) [ A , [tanh " ~ + tanh " + ^' k/ l , (4) 
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where the integration is over the Fermi surface, and kf is the Fermi momentum. Note that the kernel of the integral 
is identical to that in zero field, and that the field dependence is only in the thermal factors. Indeed, as the Doppler 
shift is the same for the quasiparticles absorbing and emitting photons the difference in energy between the initial 
and the final states remains unchanged, whereas the thermal factors, which depend on the local value of the chemical 
potential, are strongly affected by the magnetic field. 

To obtain the field dependent Raman intensity we spatially average the local response given in Eq.(^) according 
to Eq.(||). A crucial observation is that since the spatial average is performed over all directions of the superfluid 
velocity, the result does not depend on specific position ky at the Fermi surface. Consequently, spatial averaging 
decouples from the integration over the Fermi surface, and we arrive at the following universal scaling relation for the 
Raman intensity 

*"(«; H) = 0)F(Sl/2T, E H /T)/F(fl/2T, 0), (5) 



where Eh — aAo^/H/H^ is the typical quasiparticle Doppler shift, which is the energy scale introduced by the 
magnetic field, and 

F(x,y) = -[ zalz [ d6 tanh (— ~ \ f^-— cos 9) (6) 
7T Jo Jo 2 V 8 z 

is a generalization of the thermal function F(e/T,0) = tanh(e/2T). In general the function F has to be evaluated 
numerically, however it can be obtained analytically in the important limit T = when it depends only on the ratio 
x/y: 

p (v a _ / 1 - V(2^ 2 ) w>l; 
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here w = (2/Tr) 1 ^ 2 (x/y). Then the ratio of Raman intensities in any channel at low temperatures is given by 

x"t=q{^h) , n > 
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The intensity is most strongly affected for the Raman shifts below the average Doppler shift Eh, where the the 
scaling function Fq is almost linear in flH^ 1 ^ 2 , while for large frequencies the field dependent correction is small and 
linear in the applied field. In Fig.[j] we also show that the main features of the scaling function remain robust at low 
temperatures. 

So far we have made no assumptions regarding the specific symmetry of the order parameter or the particular 
Raman geometry. Assuming a 4r 2 -i/ 2 ~ wave symmetry, we obtain that the intensity at small frequencies <C Eh 
in the B\ g and Bzg channels become quartic and quadratic, respectively, compared to cubic and linear dependence 
in the absence of field. It also follows that the height of the 2Ao peak in the B^g channel decreases linearly with 
the magnetic field. Finally the integrated normalized Raman intensity df2[Ax"(fi; 0)] scales with \fTl, 

while the integral of the change in the signal itself depends on the particular Raman geometry . 

Superfluid density. Using the Green's function given in Eq.(Q) we obtain that the relative change in superfluid 
density is given by the generalized Yosida function jl3j 

5n s (T,H) 1 v d ,E k E H \ 
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where E k = -\/Ck + A 2 ., and the function F was defined in Eq.(@). It is clear from the zero-temperature limit Fq of 
the thermal function F, and from Figjl], that at low temperatures T <C Eh the energy scale determining the range 
of F(E k /T, Eh /T) as a function of E k is the magnetic energy Eh rather than T, so that the behavior of transport 
coefficients is drastically different from the H = case. As we show in Fig.^J, while in absence of magnetic field the 
superfluid density n s decreases linearly with T due to the linear low-energy density of states of a superconductor with 
lines of nodes of the gap, the low temperature behavior in the applied field becomes 

8n s (T,H) _ [8 E H , 2 V^T 2 
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and the temperature dependent term shows scaling with T 2 H~ X / 2 . In all the numerical work we have used the BCS 
value Ao = 2.14T C . The superfluid density can be extracted from either optical conductivity in the magnetic field, 
which has not been measured yet, or from the fiSR experiments, which determine the London penetration depth 
\{H) as T — > 0. The exact relationship between the penetration depth and the superfluid density in the mixed state 
is not yet clearly understood, since nonlinear [B and nonlocal |l5| effects are believed to be important, however, the 
non-linear, in H, behaviour similar to that given in Eq.([lO|) has been obtained numerically when both of these effects 
are included |16|. The inset of Fig.^| shows the data of Ref. || for a three-crystal mosaic of YBCO with our best 



linear, in \JH, fit, which gives A(0) = 1142A and H c2 /a 2 = 88T, compared to A(0) = 1155A used in Ref. || for a 
linear in H fit. Because of nonlinear effects the obtained value for H C 2 is a low-end estimate for the upper critical 
field. 

The scaling behavior given by Eq.(^|), is shown in Fig.|| where the numerical results follow the scaling curve up to 
T j 'Eh « (T / Aq)(H C 2 / H) 1 ' 2 ~ 1. For larger T/ Eh results obtained for small Eh and, consequently, low temperatures 
still scale, while for larger fields the scaling is broken since corresponding temperatures become high. 

We note that in conventional linear response theory the relative change to the spin part of the Knight shift in fields 
HbH < J 1 is also given by Eq.(^). However, in the present problem there is an additional magnetization due to the 
vortex lattice, which has to be computed from the grand potential |T^| . This contribution is small in the temperature 
dominated regime T > E H , while at low temperatures M » (a 2 ir 2 T 2 A N(0))/(6E H H c2 ). 

Spin-lattice relaxation rate. Short-range antiferromagnetic correlations, which are believed to be important for the 
Cu NMR relaxation rate in the cuprates, cancel on the oxygen sites, resulting in a normal Korringa behavior seen in 
experiment [j3||i"§|l. Then the relaxation rate T^ 1 can be calculated using the low frequency limit of the uniform spin 
susceptibility, yielding 

+ °° x d E-v s k 

iV 2 (£)— tanh * , (11) 
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where 1^ is the relaxation time at T c , and N(E) is the superconducting density of states. It is important to note that 
the NMR experiments measure the decay of magnetization M(f) oc exp(— t/Ti), rather than T\ directly, so that it is 
M(t) that has to be spatially averaged. For a distribution of 2i(r) the average magnetization cannot be described by 
a single relaxation rate, and, if such a fit is made, the obtained value of T\ is different depending on whether short or 
long time scale behavior is analyzed. At low temperatures T <§; T c , if T > Eh, the times over which measurements are 
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done, t ~ Tij correspond to short time scale in the field-dependent term, and the magnetization decay approximately 
follows a single relaxation rate behavior with 



T C T[ C) _ TT 2 T 2 nE% 
~J\T~ ~ ~3~A| + 2 Ag n 
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On the other hand, for T c 3> Eh 3> T, the single relaxation time picture breaks down completely due to strong spatial 
variations of Ti. Even though M(i) is still described by the approximate relaxation rate given by Eq.(|l2|) for t -C Ti, 
in the experimentally relevant region M(i) oc exp(-i/T 1 (H - 0) )A 1/2 - We are aware that an independent analysis of 
NMR magnetization data, including a wider range of temperatures and fields, is being carried out, using an approach 
similar in spirit to this one fujfl . 

To conclude we have presented an approach to the calculation of thermal and transport properties of d-wave 
superconductors in the mixed state over a wide range of temperatures and fields, considering the contribution of the 
extended quasiparticle states. We obtained the explicit form of the scaling functions in the low temperature regime 
T < Eh, and observed the breakdown of scaling at higher temperatures. Our results agree qualitatively with the 
fiSR measurements of the penetration depth in the vortex state. 
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FIG. 1. Zero temperature scaling function Fq(w) (solid line). The result of numerical evaluation of the function F, given in 
Eq.Q, is plotted for T = 0.03Ao for comparison (dashed line). 

FIG. 2. Superfluid density n s as a function of reduced temperature for different magnetic fields: Eh = (solid line), 
Eh = O.lTc (dashed line), 0.2T C (long-dashed line), 0.3T C (dot-dashed line). Inset: zero-temperature penetration depth data 
from /iSR experiment plotted vs. vH. Solid line: best linear fit. Slope corresponds to H C 2/a 2 = 88T. 

FIG. 3. Full numerical evaluation of the superfluid density from Eq.(^). Dashed line: low temperature result from Eq.(|lo|). 
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